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Abstract 

This work considers the g-tensor anisotropy induced by the flexural thermal vibrations in one- 
dimensional structures and its role in electron spin relaxation. In particular, the mechanism of 
spin-lattice relaxation via flexural modes is studied theoretically for localized and delocalized elec- 
tronic states in semiconducting carbon nanotubes in the presence of magnetic field. The calculation 
of one-phonon spin-flip process predicts distinctive dependencies of the relaxation rate on temper- 
ature, magnetic field and nanotube diameter. Comparison with the spin relaxation caused by 
the hyperfine interaction clearly suggests the relative efficiency of the proposed mechanism at suffi- 
ciently high temperatures. Specifically, the longitudinal spin relaxation time in the semiconducting 
carbon nanotubes is estimated to be as short as 30 fis at room temperature. 



PACS numbers: 85.35.Kt, 85.75.-d, 81.07.Vb, 81.07.Ta 



I. INTRODUCTION 



Unique electronic, mechanical and structure properties of carbon nanotubes (CNTs) con- 

n n 

tinue to be a focus of extensive investigation to date (see Refs. [Jj and|2| as well as the references 
therein). For one, it is considered to be the ultimate system for continued "scaling" beyond 
the end of the semiconductor microelectronics roadmap.- At the same time, recent studies^*^ 
drew attention to other important applications of CNTs. Their crystalline properties with 
low or no impurity incorporation allow the injection and use of electrons with polarized spin 
as an added variable for computation. Thus, CNTs are an ideal medium for the development 
of the new emerging field of spintronics.^'^ Further, the anticipated long spin relaxation times 
allows coherent manipulation of electron spin states at an elevated temperature, opening a 
significant opportunity for spin-based quantum information processing.-"^ Clearly, spin de- 
pendent properties of CNTs warrant a comprehensive investigation with the spin relaxation 
times/processes as one of the most crucial. 

Recently, electron spin relaxation in the CNTs was examined theoretically by considering 
the hyperfine interaction (HFI) with nuclear spins / = 1/2 of ^^C isotopes (with the natural 
abundance of 1.10%).— The HFI is thought to be the most important spin relaxation process 
with the role of spin-orbit coupling routinely dismissed on the ground of its weakness in the 
carbon-based structures. However, the calculation result predicts the electron spin relation 
time of the order of one second; a number much longer than that experienced experimentally. 

In the present work, we reexamine the role of spin-lattice interaction on electron spin 
relaxation in semiconducting CNTs. Although the (yf-tensor anisotropy caused by the spin- 
orbit coupling does not lead to spin relaxation, it can provide an effective mechanism when 
combined with specific thermal vibrations of the CNT (namely, the fiexural phonon modes). 
Indeed, any bending of the CNT is accompanied by local rotation of the tube axis and, thus, 
the principle axes of (7-tensor. Consequently, the electron Zeeman energy under a magnetic 
field will be affected by thermal vibrations through the local angular displacement of the 
tube axis, leading potentially to electron spin relaxation. 

The rest of this paper is organized as follows. In Sec. II, the spin-phonon interaction 
Hamiltonian is derived in terms of fiexural vibrations and the CNT (yf-tensor components. 
In Sees. HI and IV, this Hamiltonian is applied to the problem of spin-flip relaxation in a 
magnetic field for electrons localized in a finite segment of a CNT (analogous to a quantum 
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dot) and for delocalized electrons. Finally, Sec. V discusses the numerical results through 
comparison with the HFI-induced spin-flip processes. A range of parameter space is identified 
where the proposed mechanism may provide a dominant contribution. 



II. HAMILTONIAN OF THE SPIN-PHONON INTERACTION 

Consider an electron under a magnetic field B. The Zeeman Hamiltonian in an undis- 
turbed CNT can be written in the frame of reference with the axis z directed along the tube 
as 

= gnfJ'sBzSz + g±^B{BxSx + BySy), (1) 

where /ib is the Bohr magneton, S an electron spin, g± = g^x = gyy and g,^ = gzz with the 
principal values of electron gf-tensor g^x-, gyy, gzz- In the presence of fiexural thermal vibra- 
tions, the system loses the axial symmetry. If the curvature induced by fiexural deformation 
is small (compared to the CNT diameter), one can introduce local Cartesian axes x', y' , 
z' at each point along the CNT (Fig. 1). The angle between the z and z' axes is defined 
as = 9{z). Then, the Zeeman Hamiltonian in the presence of fiexural vibrations can be 
expressed in the local coordinate system (denoted as H'^) by simple substitution x, y, z ^ 
x', y', z' in Eq. ([T]) while retaining the same parameters g^^ and g^. When one transforms H'^ 
back to the laboratory coordinate x, 2, the corresponding Hamiltonian Hz shows explicit 
dependence on the angular displacement. To proceed further we fix the angle between 
X-axis and z — z' plane and take into account the smallness of Q. Keeping the linear terms 
on sin Q , the Hamiltonian can be reduced to the form Hz = h'^^ + Hs_ph with perturbation 

Hs-ph = AgiJ,B[{BxSz + BzSx) cos (p + {BySz + BzSy) sin (p\ sin 9, (2) 
where Ag = g,,- g^. 

Equation ([2]) shows that fluctuations of z' axis leads to fluctuating effective field that 
can mediate a spin relaxation. The angular displacements 6 and ip are immediate from the 
transversal deformation u, so that 

sin 6* cos V9 = — sin 6* sin (y9 = (3) 

dz dz 

In turn, the u = Vi{z) can be expressed in terms of fiexural modes. The latter are the two 
transversal acoustic modes that have orthogonal polarizations ei ± 62 (we fix their direction 



along X- and y- axes correspondingly) and possess a quadratic dispersion for mode frequency 

^x,q = Pq^ (4) 

with wave number q and polarization x = x^y}"^ Parameter (3 in Eq. (jlj) is defined by the 
diameter dt of CNT and its elastic properties, 

/3 = \dtVt, (5) 



where A = 0.56 (Refll2l) and Vt = 1.4 10® cm/s is a transverse sound velocity.— The quadratic 
dispersion law (jlj) results in singularity in phonon low energy density of states that impacts 
the dependencies on magnetic field and temperature as it will be shown below. 

If we represent the u(z) in second quantization form^^ and substitute it to the Eq. (j3]) 
and then to Eq. (j2]), after some algebra one can find 



>c=x,y q 



where gi is a linear density of CNT with length Aq, al^^ and a^^^g the operators of creation 
and annihilation of the phonon with wave number q and polarization x. 

Hereinafter we use only matrix elements of spin operators that describe the spin-flip 
transitions between the eigenstates |+) and | — ) of Hamiltonian H^^l^. The diagonal matrix 
elements can also be responsible for the phase spin relaxation providing the finite phonon 
lifetime,— we however skip this possibility assuming that phase spin relaxation is conditioned 
by spin-flip processes. 

The symmetry imposes the independence of the problem on azimuth direction angle 
of a magnetic field. For definiteness sake we fix the B -direction so that = sin a. 
By = and Bz = cos a, where a is a vectoral angle between B and the tube axis z 
(Fig. 1). Thus diagonalization of H^^^ gives rise to the |+) = (cos |, sin |) and |— ) = 
(— sin ^, cos Doing so we ignore a small difference in g,, and gj_. As a result on can find 
(- \BzSr, + Br,Sz\ +) = |5cos2a and (- {B^Sy + ByS^l +) = {Bcosa. 

The Eq. (j6j) is a basic equation that describes the electron spin interaction with flexural 
modes in CNT. It will be applied for analysis of spin-lattice relaxation in the cases of localized 
and delocalized electrons in next two sections. 
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III. SPIN RELAXATION OF LOCALIZED ELECTRONS 



Consider a localized electron in the bias potential applied perpendicular to the semicon- 
ductor CNT. Parabolic shape for such potential is a reasonable approximation if we analyze 
only the ground electronic state.-*^ That assumes quite low temperatures, ksT <^ fkvo, com- 
pared with energy space huo to the first excited state. Under these approximations the 
envelope wave function reads^ \E'o = X'^t{0'^l{v)j where C, and r] are the curvilinear coor- 
dinates associated with the circumference and longitudinal length of the deflection curve of 
CNT, X is a two-fold amount for ■j/'-function amplitudes at the A and B atoms of a primitive 
cell. We consider small bends of CNT and join rj with z axis so that longitudinal part of 
the \l/o takes the form^^ 



where 2dQ = l^f^lKfrnFuTs^ is an extension of the electron localization, m* = If? j'idt^ is an 
effective mass of the semiconducting CNT, 7 is a transfer matrix element. The implicit 
forms of the x ^-iid ipriO irrelevant to the problem under consideration. Besides, we 
ignore the modifications of the electronic states by the external magnetic field assuming that 
relevant parameter-^ {dt/2aHY is rather small {an = \Jch/eB is the magnetic length). 

The averaging of electron-phonon interaction (l6l) over the \l/o can be reduced to the 
calculation of the form-factor $(g/go) = {'^L{z)\e^''^\il)L{z)) , where go = Vs/do and 

$(a;) = e-"'. (8) 

Thus, the Hamiltonian of spin-phonon interaction Vg-ph = (^ol^l^o) for localized electron 
can be presented in canonical form 

Hs-pk = E E + (9) 

>t=x,y q 

where spin-depended operator is 



We are looking for the rate of longitudinal spin relaxation T^-^ = + w_^, which can 

be evaluated in terms of probability of spin transition from the spin states |+) 

(|— )) to the state |— ) (|+)). In turn, the w+_ and w_+ can be found by straightforward 



application of the Fermi golden rule to the Eq. (Q. The final result of such calculations 
takes the form 

27r 

M:=x,y q 

where n^^q = (aj^^a^^g) = 1/ [exp(co'^g) — 1] is the phonon population factor and huz = 
q^bB is a Zeeman splitting, a weak anisotropy of g-factor g can be neglected here. The 
evaluation of the sum in Eq. ffTTl) by means of Eq. (jlj) with matrix elements of operators in 
Eq. (fTO!) leads the final result 

where g = g,„ g±, the factor F{a) takes into account the angular dependence of longitudinal 
relaxation, 

F{a) = ^ (cos^ 2a + cos^ a) . (13) 

Equation (fT3i) describes a prominent angular dependence of spin-relaxation rate in CNT. 
It reaches the maximal value F{a) = 1 at a = 0, i.e. a magnetic field applies parallel to the 
axis of CNT, and reduces this value up to F{a) = 0.22 at a = 52.5° (Fig. 2). In the case of 
perpendicular orientation, a = 90°, only modes with x polarization contribute to relaxation, 
i.e. F{a) = 0.5. 

IV. SPIN RELAXATION OF DELOCALIZED ELECTRONS 

Let us consider the spin-fiip relaxation for delocalized electrons in semiconductor CNT. 
This case is most relevant to the spintronic applications, so we extend the consideration 
up to room temperatures. Assuming small displacements u we describe the longitudinal 
component of orbital electronic states in the vicinity of the K valley in terms of plane waves 
with wave number k counted off from the K point of the Brillouin zone, = e^^^ / \/A^. 

This approximation immediately results in momentum conservation when interaction with 
phonons describes the Hamiltonian ([6]), i.e. (fc'|e*^^|fc) = {271 / AQ)5{q — k' + k) and 

{k'\Hs-,H\k) = J2 Vk'ki^^,-^k + aL,Afc), (14) 

M:=x,y 
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where Ak = k' — k, and 



=^sign(AA:)^^^-A_(5^^^ + B^5.). (15) 

The electron scattering in the vicinity of the K' valley gives rise to the same results while 
intervalley scattering is negligible because it entails the excitation of phonons with giant q. 

The probability w'^^ of spin-flip |+) |— ) due to electron scattering k ^ k' on the 
phonons can be found in terms of Fermi golden rule. The averaging over the phonon thermal 
distribution results in 



C [n^k5 {AE -hujz- hf3Ak^) + {uAk + 1) 5 {AE - hujz + hf3Ak^)] (16) 



where C = iTF{a){AgfiBBy /2giAoP, AE = h^{k''^ -k'^)/2m* . Similarly one can obtain the 
probability of spin flop w^^. The net result for electron spin relaxation ?iJ+_ assumes the 
averaging of the w^'_ over the electron initial state occupation numbers and summation 
over non-occupied final states k' \ —), 



w+_ = 

k,k 



The summation over k, k' has to be replaced by integration in the common fashion. In the 
case of non-degenerate electrons, 



One can see that the Eqs (fT6l) - f|T8|) are conveniently expressed in terms of new variables 
u = k + k' and v = k — k' so that arguments of 5-functions in Eq. fll6p transform to the 
linear form on u. The integration over u in Eq. ( ITTl) becomes trivial and gives rise to final 
expression for the spin relaxation rate T^^ = + in the form 

= 4^ f ^) ' ^^^pk' CO* I cosh pih), (19) 



57r 

where we introduce the dimensionless parameters h = Hwz/kBT and e = 2m*P/h, and an 
integral function 

m = exp --—x - h'- -— — — . (20) 



Ae Ax J a;sinh(a;/2) 
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Taking into account Eq. ([5]) and definition of m*, parameter e can be estimated as a number 
e = 0.02 independent on dt. 

The integrant in Eq. (1201) appears as a result of the change of variable vhj x = h(3v'^ / ksT. 
Asymptotic behavior of the Eq. fl2Ul) is I{h ^ 0) ^ S/eh"^, i.e. the integral diverges at zero 
magnetic field {h 0). Actually this singularity has not physical meaning. If electronic 
scattering is taking into account, the wave numbers k and k' are not longer accurate quantum 
numbers and 5-functions in the Eq. f|T6|) transform to delta-like functions with finite width 
and maximum. This limits the value of integrant in Eq. (!20|) and function I{h) at small x and 
h correspondingly. Similar situation was discussed in Refs^*^. Practically, the singularity 
treatment means that once Zeeman energy hwz becomes of the order of electron energy 
broadening F, the integral will reach a saturation, i.e. if 

< r, (21) 

then I (huJz/kBT) I (T/kBT). As a result, zero magnetic field quenches spin relaxation 
mechanism under consideration, T^^{ljz ^ 0) ~ [fujOz)^ I (y^f^ ~^ 0- We will not consider 
this effect in more detail. Fig. 3 presents the function h^I{h) calculated with Eq. (120|) 
without considering the effect of finite F. Note that the h'^I{h) exponentially decreases in 
the limit of large h. 

V. NUMERICAL EVALUATION AND DISCUSSION 

Equations ( |T2l) and ( |T9l) define the spin relaxation rate via flexural phonon modes for 
localized and delocalized electrons. Alternative approach based on the electron spin scat- 
tering on disordered nuclei spins / = | of ^^C isotopes due to hyperfine interaction was 
done in the RefJlOi. Comparing these mechanisms is performed for chirality vector (8,4) 
(respectively CNT diameter and effective mass are dt = 0.839 nm and m* = 0.186 mo, ttlq 
the free electron mass)^° with following parameters: Ag = 0.005, g = 2.015,^^ gi = TTdtQ2 
where Q2 = 9.66 x 10^^ g/cm^ is the mass of graphene sheet per unit area.- Hereinafter all 
calculations are performed at B directed along the tube axis, i.e. a = 0. 

Figure 4 illustrates the temperature and magnetic field dependencies of spin relaxation 
rate of the electron localized at 2dQ = 30 nm. The temperature dependence is typical for 
one-phonon process of spin- lattice relaxation in semiconductor quantum dot.^^ At relatively 
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high temperatures [ksT > tujJz) the population factor for resonant phonons describes the 
direct process ~ T while the temperature independent phonon radiation dominates at 
ksT < hwz- Magnetic field dependence runs from zero through the maximum, which is 
caused by form-factor at relatively short wavelength of resonant phonon. Dependence on 
CNT diameter (Fig. 5) stems from the explicit dependencies of form-factor, the fi [Eq. 
and pi on dt- Both Figs. 4 and 5 display rather long (around few seconds) spin relaxation 
of localized electrons in CNT at low temperature that makes it attractive for quantum 
computing application. 

The magnetic field and temperature dependencies of spin-relaxation rate calculated with 
Eq. f ll9p is shown in Fig. 6 for the delocalized electron (i.e., without taking into consider- 
ation the electron scattering effects). Some decreasing attributes to only low temperature 
relaxation in actual range of the magnetic fields while no significant dependence on B re- 
veals the Tj~^ starting with T = 40 K. For the room temperature, Eq. (fT9l) predicts the 
relaxation time about Ti ~ 30 /xs in the wide range of the magnetic fields. The tends to 
increase as the temperature goes up while mechanism due to hyperfine interaction possesses 
the opposite tendency. The increase of the tube diameter suppresses spin relaxation in the 
case of both mechanisms (Fig. 7). Note that with the exception of very low magnetic field, 
which corresponds to Eq. (!2T]) . the spin relaxation via fiexural modes dominates over the 
processes of electron scattering at nuclei spins. This deduction is based on the comparison 
of present results and theory developed in Ref. 
specified, ahf/27fh= -35.8 MHz.i^ 



10, where the constant of HFI has been 



VI. CONCLUSIONS 



A mechanism of electron spin relaxation in carbon nanotubes caused by anisotropy of 
g-tensor and fiexural phonon modes is studied. Relaxation time for localized electron is 
estimated around a second at low temperature while the spin relaxation of delocalized elec- 
trons can reach few tens microseconds at room temperature. As it turned out the proposed 
mechanism is essentially more efficient than the processes of spin scattering due to hyperfine 
interaction with isotopes ^^C. This holds true up to very weak magnetic fields when Zeeman 
splitting comes up with electron energy broadening. The mechanism reveals specific depen- 
dencies on magnetic field orientation (common for both localized and delocalized electrons) 
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and on the magnetic field strength (with perceptible maximum for localized electron and 
almost fiat dependence for delocahzed one). Besides, in contrast to processes of scattering 
on nuclei spins, the temperature augmentation sufficiently increases relaxations rate. These 
particularities will facilitate the experimental recognition of the mechanism of spin-lattice 
relaxation in CNT. 
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CNT 



FIG. 1: The position of the laboratory coordinate system x, y and z, the local coordinate x', y' 
and z' related to CNT flexure and the angles between axes. The direction of a magnetic field B 
determines the angle a. 
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FIG. 2: Function graph of F{a), which defines the dependence of spin relaxation rate on the angle 
a between the direction of magnetic field and the axis of CNT for both localized and delocalized 
electrons. 
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Temperature (K) Magnetic field B (T) 

FIG. 4: Calculated spin relaxation rates of localized electron in a (8,4) CNT as a function of 
magnetic field and temperature. Curves 1, 2, and 3 are obtained for B = IT, 0.3T, and O.OOIT 
correspondingly. Curves 4 and 5 are plotted for T = lOK and IK. The size of electron localization 
is 2do = 30 nm. 
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FIG. 5: Calculated spin relaxation rates of a localized electron as a function of CNT diameter 
for different values of magnetic field and temperature and temperatures T = IK (left panel)and 
T = 10 K (right panel). Curves 1, 2, and 3 correspond to B = O.OOIT, IT, and 3T respectively. 
The size of electron localization is 2do = 30 nm. 
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FIG. 6: Calculated spin relaxation rates for mechanism caused by flexural modes (solid lines) and 
due to electron scattering on ^^C isotopes (dashed line) for delocalized electron in a (8,4) CNT as 
a function of magnetic field at T = 4K (left panel) and temperature (right panel). The curves 2, 
3 and 4 are obtained for B = 0.05T, IT, and 5T correspondingly. Solid line 1 practically does not 
change in range of the magnetic field from B = 0.05T up to B = 5T. 
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FIG. 7: Calculated spin relaxation rates of delocalized electron as a function of CNT diameter for 
different values of magnetic field and temperature. Contributions of both HFI (dashed line) and 
spin-phonon interaction (solid line) are presented for the purpose of comparison. Curves 1, 3, and 
6 correspond to B = 0.05T. Curves 2, 4, and 7 correspond to B = 2T. Curve 5 does not change in 
the range of the magnetic field B = 0.05 — 2T. 
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